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Abstract 

Let Vr be a lattice periodic potential and A and (j> external electromag- 
netic potentials which vary slowly on the scale set by the lattice spacing. 
It is shown that the Wigner function of a solution of the Schrodinger equa- 
tion with Hamiltonian operator H = ^(— iVx — A{ex))^ + Vr{x) + (f>{£x) 
propagates along the flow of the semiclassical model of solid states physics 
up an error of order e. If e-dependent corrections to the flow are taken 
into account, the error is improved to order e^. We also discuss the prop- 
agation of the Wigner measure. The results are obtained as corollaries of 
an Egorov type theorem proved in IPST3I . 

1 Introduction 

One of the central questions of solid state physics is to understand the motion 
of electrons in the periodic potential which is generated by the ionic cores. 
While this problem is quantum mechanical, many electronic properties of solids 
can be understood already in the semiclassical approximation |AsMel IKol IZa| . 
One argues that for suitable wave packets, which are spread over many lattice 
spacings, the main effect of a periodic potential Vr on the electron dynamics 
corresponds to changing the dispersion relation from the free kinetic energy 
Efieeip) = ^P^ to the modified kinetic energy En{p) given by the n**^ Bloch 
function. Otherwise the electron responds to slowly varying external potentials 
A, (f> as in the case of a vanishing periodic potential. Thus the semiclassical 
equations of motion are 

r = \/En{K), k = -V(t){r) +r X B{r) , (1) 

where n — k — A(r) is the kinetic momentum and B = curM is the magnetic 
field. (We choose units in which the Planck constant h, the speed c of light, and 
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the mass m of the electron are equal to one, and absorb the charge e into the 
potentials.) The corresponding equations of motion for the canonical variables 
(r, k) are generated by the Hamiltonian 

H,,{r,k)^E^{k~ A{r))+4>{r), (2) 

where r is the position and k the quasi-momentum of the electron. Note that 
there is a semiclassical evolution for each Bloch band separately. The distinction 
between the canonical variable fc, the Bloch- or quasi-momentum, and the ki- 
netic momentum k = k—A(r) is often not made explicit in the physics literature. 
It is, however, crucial for the formulation of the precise connection between the 
semiclassical equations of motion ^ and the underlying Schrodinger equation 
©■ 

we use adiabatic perturbation theory in order to understand on a 
mathematical level how these semiclassical equations emerge from the underly- 
ing Schrodinger equation 

15, ^(y, s) = (i( - iV, - A{ey)f + Vriy) + ^(ey)) V(y, s) (3) 

in the limit e ^ at leading order. In addition, the order e correction to ^ 
are established, see Equation {Tj). 

In (PJ the potential : R'' — > M is periodic with respect to some regular 
lattice r generated through the basis {71, . . . , 7^}, 7^- S M'', i.e. 

r = |a; e M'* : a; = ^^=1"^ lo some a G Z^'j 

and Vr( • + 7) = Vr(-) for all 7 G F. The lattice spacing defines the microscopic 
spatial scale. The external potentials A{ey) and (piey), with A : R'^ ^ R'^ and 
(/) : R'^ — > R, are slowly varying on the scale of the lattice, as expressed through 
the dimensionless scale parameter e, e ^ 1. In particular, this means that the 
external fields are weak compared to the fields generated by the ionic cores, 
a condition which is satisfied for real metals even for the strongest external 
electrostatic fields available and for a wide range of magnetic fields, cf. |AsMe| . 
Chapter 12. 

Note that the external forces due to A and (j) are of order e and therefore 
have to act over a time of order e^^ to produce finite changes, which is taken as 
the definition of the macroscopic time scale. Hence, one is interested in solutions 
of for macroscopic times. The macroscopic space-time scale (x, t) is defined 
through X = ey and t = es. With this change of variables Equation (j2Jl reads 

iedtnx,t) = (^^{-isV.,-A{x))^ + Vr{x/e) + q}ix)^iP%x,t) (4) 

with initial conditions ^p^{x) = e~'^^'^^p{x/e). If Vr — 0, then the limit e ^ in 
Equation |0J is the usual semiclassical limit with e replacing h. 

The problem of deriving |^ from the Schrodinger equation Q in the limit 
e ^ has been attacked along several routes. In the physics literature is 
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usually accounted for by constructing suitable semiclassical wave packets. We 
refer to |Lul IKollZaj . The few mathematical approaches to the time-dependent 
problem extend techniques from semiclassical analysis, as the Gaussian beam 
construction [UnTlinUH] . or Wigner measures GMMP, BFPR, BMP . 

In this note we explain and elaborate on recent results from jPSTaj. In 
IPST3I we derived ^ from for quite general external potentials A and (j). 
The construction is based on the space-adiabatic perturbation theory developed 
in IPSI'iin^ . see also |MeSo| and the contribution of G. Ncnciu in the present 
volume. The crucial observation is that the step from © to (^l involves actually 
two approximations. Semiclassical behavior can only emerge if a Bloch band is 
separated by a gap from the other bands and thus the corresponding subspace 
decouples adiabatically from its orthogonal complement. The dynamics inside 
this adiabatic subspace is governed by an effective Hamiltonian h^g, which is 
explicitly given as an e-pseudodifferential operator. Eventually, the semiclassical 
limit of hlff leads to (0). 

Hence ^ needs to be reformulated as a space-adiabatic problem. This has 
been done first in (HST, for the case of zero magnetic field and then in 1PST3I for 
general electric and magnetic fields. The results obtained in this way constitute 
not only the derivation of the semiclassical model ^ in this generality, but they 
allow to compute systematically higher order corrections in the small parameter 
e. It turns out that the electron acquires a fc-dependent electric moment An{k) 
and magnetic moment A4„(/c). If the -n}^ band is nondegenerate (and isolated) 
with Bloch eigenfunctions ipn{k,x), the electric dipole moment is given by the 
Berry connection 

Anik) =i{Mk),^Mk)) ■ (5) 
and the magnetic moment by the Rammal- Wilkinson phase 

Mnik) = ^{VMk), x(iJpe,(fc)-£;(fc))V7A„(A:)). (6) 

Here (• , •) is the inner product in L'^{R'^/T) and Hpcr{k) is of © with <f) = 
= A for fixed Bloch momentum k. Note that En, An and Ain are r*-periodic 
functions of k, where F* is the lattice dual to F. Hence one can as well think of 
them as functions on the domain M* = R'^/F*, the first Brillouin zone. 

The semiclassical equations of motion including first order corrections read 

f ^ \7,,(^En{K)-eB{r)-Mnin)j -ekx^nin), 

(7) 

k = -Vr(^(l>{r) - eB{r) ■ Mnin)) +r X B{r) . 

with ri„(fc) = V X An{,k) the curvature of the Berry connection. 

In order to state the precise connection between the semiclassical equations 
of motion resp. their refined version lO and the underlying Schrodinger 
equation Q), we need some more notation. Let 

i/^ = i( - ieV, - Aix))^ + Vrix/e) + ^{x) (8) 
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be the Hamiltonian of Under the following assumption on the potentials, 
which will be imposed throughout, is self-adjoint on H'^iW'-). Here C^(M'^) 
denotes the space of smooth functions which are bounded together with all their 
derivatives. 

Assumption. Let Vr he infinitesimally bounded with respect to —A and assume 
that (f) e C^{W\M.) and Aj G C^(R'^,R) for any j G {1, . . . ,d}. 

To each isolated Bloch band En there corresponds an associated almost 
invariant band-subspace n^L^(R'^). The orthogonal projector onto this 
subspace is constructed in IPST3I . Only for states which start in this subspaces 
and thus, by construction, remain there up to small errors, the semiclassical 
equations of motion ((TJ can have any significance. 

The flow of the dynamical system (Q is denoted by $* 
canonical coordinates (r, k) — [r, k + A(r)) by 



T>2d 



T>2d 



The existence of the smooth family of diffeomorphisms is not completely 
obvious from {T)) alone, but follows from the Hamiltonian formulation of iQ 
presented in the next section. 

Notation. Throughout this paper we will use the Frechet space 

C = C{J°(R2''), 

equipped with the metric dc induced by the standard family of semi-norms 

l|a|U = l|5"a|U, aeNl'', 
and the subspace of r*-periodic observables 



per 



{a G C : a(r, fc + 7*) = a{r, k) V7* G T*} 



We abbreviate dc{a) := dc(a,0). 



The main result of IPST3I on the semiclassical limit of (@J is the following 
Egorov type theorem. 



Theorem 1. Let he an isolated, non- degenerate Bloch hand. For each finite 
time-interval / C M there is a constant C < 00, such that for all a G Cpor with 
Weyl quantization a = a{x, — ieV^:) one has 



^iHH/e^^-iH't/e 



a o $ 



H^ 



< eCdc(a) 



(9) 



B{L' 



and 



<e^Cdcia). 



(10) 



e(L2(Rd)) 
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Remark. The corresponding statement in 'PST3' does not make explicit the 
dependence of the error on the observable a. However, the more precise version 
formulated here is a standard consequence of the Calderon-Vaillancourt theorem 
and the fact that composition with is a continuous map from C into itself. <)> 

Remark. On an abstract level the distinction between the functions <i>* and 
is immaterial, since both functions express the same dynamical flow in two 
systems of coordinates. However, the distinction between the systems of co- 
ordinates becomes important when the quantization is considered. The Weyl 
quantization appearing in and Hl()(l must be understood with respect to the 
system of coordinates (r, k) . Analogous consideration hold true for formulas 
involving a Wigner transform, as in Corollary [21 (} 

The main objective of this note is to elaborate on Theorem ^ in order to 
make contact to alternative approaches and results on the semiclassical limit of 
(0J. This are, as mentioned above, Wigner functions '( IMMl^l IEeWI IEMF| . 
semiclassical wave packets jLul IKol IZal ISuNi) and WKB-type solutions of (0} 
| Bul CRT, DGJH. We focus on the semiclassical transport of Wigner functions 
and Wigner measures in the following. Before we do so, it is worthwhile to first 
examine the equations of motion jT)) in some more detail. 

2 The refined semiclassical equations of motion 

The dynamical equations l(7)), which define the e-corrected semiclassical model, 
can be written as 

r = \7 K.Hscir, k) - e k X Qn{K) , 

(11) 

k = —\/rHsc{r,K)+rxB{r) 

with 

i7sc(r,K) ■.^Er,{K) + (b{r)~eMn{K,)-B{r). (12) 

We shall show that are the Hamiltonian equations of motion for H12|) with 
respect to a suitable e-dependent symplectic form Qb.c- The semiclassical equa- 
tions of motion ijJJ) are defined for arbitrary dimension d. However, to simplify 
presentation, we use a notation motivated by the vector product and the duality 
between 1-forms and 2-forms for d = 3, which we briefly explain. 

Notation. If d 7^ 3, then B, f7„ and A^„ are 2-forms with components 
By (r) = d,Aj (r) - d,A, (r) , 17,, (fc) = d,Aj (k) - djA^ik) 

and 

M^,{k) = Re i (a,^„(fc), (i/p„. - E){k) d,Mk)) ■ 

For d = 3 a, 2-form Bij{r) is naturally associated with the vector Bk{r) = 
£kij Bij (r) . We use the convention that summation over repeated indices is 
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implicit. Then in iQ the inner product B ■ Mn refers to the product of the 
associated vectors and we generaUze the notation to arbitrary dimension d using 
the inner product of 2-forms defined through 

d d 

B-M:= *-^{B A *M) ^Y^^BijM^j , 

3 = 1 »=1 

where * denotes the Hodge duahty induced by the euchdian metric. In the same 
spirit for a vector field w and a 2-form F the generalized "vector product" is 

d 

{w X F)j := {*-^{w A *F))j = ^w,Fij , 

1=1 

where the duality between 1-forms and vector fields is used implicitly. <^ 

We keep fixed the system of coordinates z = (r, k) in R^** for the following. 
The standard symplectic form 80 — Oo(^)/m dzmAdz;, where l^m G {1, . . . , 2d}, 
has coefficients given by the constant matrix 

where I is the identity matrix in Mat((i, R). The symplectic form, which turns 
into Hamilton's equation of motion for H^c, is given by the 2-form Qb,s = 
QB,e{z)im dzm A dzi with Coefficients 

For £ = the 2-form 0_B,e coincides with the magnetic symplectic form Qb 
usually employed to describe in a gauge-invariant way the motion of a particle 
in a magnetic field f |MaRa| . Section 6.6). For e small enough, the matrix H13|l 
defines a symplectic form, i.e. a closed non-degenerate 2-form. 

With these definitions the corresponding Hamiltonian equations are 

eB,e{z) Z^dHsciz), 

or equivalently 

f B{r) -I \fr\_f VrH{r, k) \ 
\ I e r!„(«) J\k) y y.H{r,K) J ' 

which agrees with . We notice that this discussion remains valid if ri„ admits 
a potential only locally, as it happens generically for magnetic Bloch bands. 

The symplectic structure is therefore determined by the magnetic field B{r) 
and by the curvature of the Berry connection f2(A;), which encodes relevant 
information about the geometry of the Bloch bundle ipnik, ■) ^ k G M* . One 
can show that, whenever the Hamiltonian has time-reversal symmetry one 
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has that ^ln{—k) = — r2„(fc). Moreover, if the lattice F has a center of inversion, 
then ri„(— fc) = r2„(fc). Thus, the two symmetries together imply that fln{k) 
vanishes pointwise. But there are many crystals which do not have a center of 
inversion and, more important, in the presence of a strong uniform magnetic 
field the time-reversal symmetry is broken. The latter is the typical setup to 
describe the Quantum Hall Effect, a situation in which the curvature of the 
Berry connection plays a prominent role. Indeed, the equations of motion lO 
provide a simple semiclassical explanation of the Quantum Hall Effect. Let us 
specialize {Tj) to two dimensions and take B{r) — 0, (p{r) = —£ ■ r, i.e. a weak 
driving electric field and a strong uniform magnetic field with rational flux. 
Then, since k — k, the equations of motion become r = VkEn{k) + £-^^ln{k), 
k — £, where i7„ is now scalar, and £^ is £ rotated by tt/2. We assume initially 
fc(0) — k and a completely filled band, which means to integrate with respect 
to k over the first Brillouin zone M*. Then the average current for band n is 
given by 

Jn = I dkr{k) = f dk (Vfc£;„(fc) - £^nn{k)) = -£^ f dkn^ik) . 

dk^nik) is the Chern number of the magnetic Bloch bundle and as such 
an integer, cf. |TKNN| . Further applications related to the semiclassical first 
order corrections are the anomalous Hall effect (JNM, and the thermodynamics 
of the Hofstadter model |GaAv| . 

3 Semiclassical transport of Wigner functions 

Theorem n provides a semiclassical description of the evolution of observables. 
The most direct way to turn it into a description for the semiclassical evolution 
of states is via duality, i.e. via the Wigner function. Recall that according to 
the Calderon-Vaillancourt theorem there is a constant C < oo depending only 
on the dimension d such that for a G C one has 

KV', aV)L2(R-) I < Cdc{a)M\ (14) 

Hence, the map C 9 a t-* (V', aip) G C is continuous and thus defines an element 
wf of the dual space C', the Wigner function of tp- Writing 

{■ip,a^)=:{wf,a)c',c=- dqdp a{q,p) wf {q,p) (15) 

and inserting into IjlSI) the definition of the Weyl quantization for a G iS(M^'') 

(a^){x) = -l- f dedya(i(x + y),ee)e'«-(^-«)V(y), 

one arrives at the formula 

^t{<l.P)^7ihz I dCe'«-PV*('? + e^/2)V'('?-£C/2) (16) 
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for the Wigncr function. Direct computation yields 



|w^||L2(K2d) =e 



{2n) 



-d/2 



Therefore, wf G L^(R^'') for aU e > 0, which explains the notion of Wigner 
function. Although wf is obviously real-valued, it attains also negative values 
in general. Hence, it does not define a probability distribution on phase space. 
However, it correctly produces quantum mechanical distributions via H15|) . 

With this preparations we obtain the following corollary of Theorem^ which 
says that the Wigner function of the solution of the Schrodinger equation is 
approximately transported along the classical flow of ^ resp. Q). 

Corollary 2. Let En be an isolated, non-degenerate Block band. Then for each 
finite time-interval / C K there is a constant C < oo such that for t E I , a E Cpcr 
and for V'o G HfjL^(R'') one has 



and 



Here ipt 



(( 



,,'00 



C',C 



C',C 



<sCdc{a)\\M' 



<e^Cdcia)Uo\ 



t/e^p ig iliQ solution of the Schrodinger equation 

Remark. When proving results for the transport of Wigner functions or Wigner 
measures it is common, e.g. |GMMPI IMMP BMP|. to write down the trans- 
port equation for We{t) := wf'^ o $q * instead of using the flow <i>g. Clearly 
our results can be reformulated in this way, cf. Corollary but the resulting 
transport equation looks complicated compared to the simple dynamical system 
Q governing its characteristics. <^ 

Proof of Corollary\^ The result is rather a reformulation of Theorem ^ than a 
real corollary. According to the Definition 1151) and Theorem Q] one has 



' ae il„'0o)L2(Ri) 



= (V'o, ao$*?/'o)L2(Hd) -|-0(e2). 

Since the map C 3 a t-^ a o eCis continuous, the duality relation 1)1 5|l can 
be applied again and yields 



□ 
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Since the functions En, Ain and r2„ appearing in the equations of motion |(7J) 
are all F* periodic, the natural phase space for the flow lO is X T* rather than 
Here T'' /T* is the first Brillouin zone M* equipped with periodic 

boundary conditions. Hence one can fold the Wigner transform onto the first 
Brillouin zone and define 



w:^^^{r,k)^ ^ u;f(r,fc + 7*) for (r, fc) G R"^ x T'' . (f7) 



e red V 



Then for periodic observables a it follows that 

/ drdka{r,k)wf{r,k) = I dr dk a{r, k + wf (r, k + ^*) 

= ^ J dr dk a{r, k) wf {r, k + J*) 

^ /drd*„(r..,<,(r,.)^ 

Thus the statement of Corollary |21 in terms of the reduced Wigner function 
becomes 

(^t, aiJt)mR^) - / dr dfc a(r, k) (wf^^^ o $7*) (r, k) + 0{e') . 

Note that the reduced Wigner function wfred coincides with the "band- Wigner 
function" of [MMP, and the "Wigner series" of [BMP], both defined as 

This follows by a simple computation on the dense set V' G iS(M''): 

E / dee'«-^*e'«-'=^(r + ee/2)^*(r-ee/2) 



\M* 



= ^E^''''^('^ + ^^/2)'/^*('^-e7/2), 
where (Sr(C) = S7er ^ t)- We used the Poisson formula 



9 



4 Classical transport of the Wigner measure 

We now turn to the Wigner measure. Recall that the Wigner function ■wf{q,p) 
can be negative and, as a consequence, does not define a probability distribution 
on phase space. In the limit e —> however, wf weakly converges to a positive 
finite Radon measure fi^ e A^^(R^'') on phase space M^'^, the Wigner measure 
of V'- For surveys on Wigner measures see e.g. [LiPal IUMMPI . 

Proposition 3. Let Sj ?° and {ipj}jeN C L'^{W^) be bounded, then the set 
{wfj }jGN C C is weak-* compact and every limit point pL ^ C defines a bounded 
positive Radon measure, called a Wigner measure o/jV'jljeN- 

Proof. The Calderon-Vaillancourt theorem (|14|) implies that {w^/} C C is 
bounded. Hence, it is weak-* compact. By (|15|) and the semiclassical sharp 
Garding inequality, e.g. Theorem 7.12 in |DiSj| , it follows that for each a > 
there is some C < oo such that 

{wt, a)c'fi >-Ce for aU?A e L^{R'^) . 

This implies the positivity of all limit points in C, which therefore define mea- 
sures. 

Let /i g C be such a limit point with, after possible extraction of a subse- 
quence, wfj ^ fi. From (|15() it follows that 

{wf, l)c'fi = MUm^) for aU^ e L^R") , 

and thus, 

^{R^'^) =^ sup{^{K) ■.KcR^'^ compact} 

< (m, 1)c',c = lim {wfj, l)c',c = lim HV'j Ilis.Rd) . 

Hence, /x is bounded. □ 

However, not all limit points are physically sensible. For example, the 
bounded sequence tpj{x) := ^o{x — j) G L^(R) has a limit point in C, some 
Banach-limit type functional, but the corresponding measure is zero. More gen- 
erally, there are many continuous linear functionals on C which are zero on the 
(non dense) subset C^{R'^''^). 

Definition. A sequence {ipj}jeN remains localized in phase space (with respect 
to {sjljeN), if it is compact at infinity, i.e. 

lim lim sup / dx \^jjj{x)\'^ — 0, 

j^oo J\a:\>n 

and e-oscillatory, i.e. 



lim lim sup— / dp \'il!j{p/ej)\'^ ~ 0. 

n^oo j^oo £j J\p\>n 
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Proposition 4. Let wf- fJ. in C with {V'jljeN C L'^{W^) bounded and local- 
ized in phase space, then /i has total mass 

^(R2<i)^ lim (18) 

and its marginals are given through the weak limits (in -Ad^ ) of the quantum 
mechanical distributions, i.e. for all a G CfKM"^) one has 

^i{dq,dp)aiq) - lim [ dq\^,iq)f aiq) , (19) 

H{dq,dp)a{p) = lim ej'^ / dp |V'j(p/ej)P a(p) . 



Proof. We start with the position marginal (|19|l . Let a G C|if'(R'') and let 
{an}«eN C C^{M.'^) and {XnjneN C C^(M'') satisfy a^iq) = a{q) and Xn{p) = 1 
for |g| < n resp. \p\ < n. Then, by dominated convergence, 

fi{dq,dp)a{q) = lim / fi{dq,dp) an{q)Xn{p) 

= lim {fi, anXn)c'.c = 1™ lim {wf', anXn)c',c 

n — ^oo ■ n — *oo J — >oc -' 



= lim / dq\il;j{q)\ a{q) + R , 

]^oo J 

where 

\R\ < lim lim \{wf^ , {a - anX7i))c'.c\ 

n^oo j — >oo 

< hm hm ( | (t/'j , {a ~ aoOn) ijjj ) | + | (V'j , {aoOn ~ a^Xj,) 

= lim lim (KV'j, (a -ax^) + {{ipj, (ax^ -a^X^)V'j)l) 

< lim lim (||a-!/;j|| 11(1 - Xn)V'j|| + - an)V'ill llXnV'jII) 
= 0. 

For the last equality we used that {ipj} is localized in phase space. In order 
to prove ((T^ also for a G C° note that we just proved that the right hand 
side of H19|l defines a measure. Hence, the result follows again by dominated 
convergence. The statements about the momentum marginal and the total mass 
follow analogously. □ 

We now turn to the propagation of Wigner measures. As remarked in the 
introduction, a popular approach to the semiclassical limit of Q is to determine 
the resulting transport equation for the Wigner measure associated with an e- 
dependent initial condition 

Corollary 5. Let En be an isolated, non-degenerate Bloch hand. Let /iq be 
the Wigner measure of a bounded sequence {"00,^} with i/jqj G Iln L'^{R'^), i.e. 
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Then the Wigner function wf^'^ of the the time- evolved sequence 

has the weak-* limit G C'p^^. given through 

/it = A^o o ^o"* ■ (20) 
In particular, /it is a positive bounded measure and solves the transport equation 

fi + VEn{k - A{r)) ■ V,/i - (V(/.(r) - diEn{k - A(r)) Vv4,(r)) • Vfe/t = 

in the distributional sense. 

Similar results were proved in |MMPI IGMMPI |BFPR" for the case of van- 
ishing external potentials A and (j)- For vanishing magnetic potential A but 
nonzero electric potential (f> they follow from the results in |HST| or ^MP^. 

Proof of Corollary\^ According to Corollary [21 we have for a G Cpor that 



((<---<o.^oci>r),a)^,^ 



<SjCdc{a)UoJ^- 



Taking the limit j ^ oo on both sides yields the existence of the limit fit and 
at the same time H2()|l . The transport equation for /it follows by taking a time- 
derivative in (|2U|) and recalling that $o is the Hamiltonian flow of ([SJl. □ 
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